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CHAPTER I 


INTRODUCTION 

During a radar cross section (RCS) measurement, the target is held in place 
using a support structure which can affect the accuracy of the measurement, es- 
pecially when the target under test has a relatively low RCS. Therefore, a support 
structure which has a low RCS is desired for scattering measurements. Among 
the commonly used structures are metal ogival pedestals, styrofoam columns and 
dielectric straps. These three structures are shown in Figures 1 through 3, respec- 
tively. 

The metal ogival pedestal is commonly used as a target support because it can 
hold heavy targets and has a relatively low RCS, especially at high frequencies. 
Lai and Burnside [1,2] have shown that the RCS of a metal ogival pedestal is 
directly proportional to the square of the wavelength of the incident field and is 
approximately 20 dB below a square meter at 1 GHZ. Thus, metal ogival pedestals 
do not appear to be a good choice for supporting targets at low frequencies; i.e. 
below 1GHz. For a horizontally polarized incident wavefront the RCS of styrofoam 
columns is approximately equal to that of metal ogival pedestals while for vertical 
polarization the RCS of styrofoam columns is higher [2]. Furthermore, to support 
heavy targets the volume of the styrofoam columns should be rather large. Thus, 
styrofoam columns do not appear to be a good choice either. Gupta, Lai and 
Burnside [2,3] studied the scattered fields from dielectric straps and concluded that 
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Figure 2: (a)Styrofoam column. (b)Target on a styrofoam column. 
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Figure 3: (a) A Dielectric strap, (b) Target attached to straps. (c)Target hung 

by straps. 


they perform much better at lower frequencies [2,3]. One advantage of the straps is 
that when the plane wave illuminates the straps at or near edge-on incidence, the 
RCS of the strap is quite low at low frequencies. Even though its RCS increases 
when the frequency increases, it is still lower than that for metal ogival pedestals. 
As a result, it is possible to align the straps with the low RCS direction of the 
target, so the minimum RCS of the target coincides with the minimum RCS of the 
strap. In addition, the straps can be built using composite material so they can 
handle very heavy targets. Therefore, the dielectric straps appear to be a good 
choice for supporting targets at low frequencies. 

The straps can be used to hold a target using at least two basic methods. 
The first one is by connecting the straps directly to the target as shown in Figure 
3(b). The second method is by wrapping the strap around the target as shown 
in Figure 3(c). The second method appears to be easier in that the target does 
not have to be modified. Both methods are expected to have a relatively low RCS 
when used as a target support. If the strap is directly attached to the target, 
the RCS of the support will be dependent on how well the strap is connected 
using a low scattering attachment. One the other hand, the second approach 
is dependent on the scattering level of the strap as it surrounds the structure. 
In order to evaluate the scattering level for this type of attachment, this report 
will determine the scattering performance of a thin dielectric strap surrounding a 
perfectly conducting structure. 

The method of analysis can be outlined in the following steps. First, the cur- 
rent induced within the strap by an incident plane wave is found by employing the 
method of moments and the Uniform Geometrical Theory of Diffraction (UTD). 
After the current is found, the UTD is used in writing an integral for the radiated 
field based on the strap currents. The integral is simplified by using a station- 
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ary phase point integration method. Then, the resulting expression is evaluated 
numerically. This method is used to evaluate the RCS of a strap surrounding a 
perfectly conducting structure, and reasonable agreement was found between the 
computed and the measured results. 

The different chapters of this report are as follows. In Chapter II, a method is 
developed to compute the scattered fields from a thin dielectric strap surrounding 
a perfectly conducting convex surface. In Chapter III, the method is employed to 
solve for the backscattered field from a strap surrounding a perfectly conducting 
convex structure. Then, in Chapter IV some computed as well as measured results 
are presented. The conclusions are given in Chapter V. 
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CHAPTER II 


SCATTERING FROM THIN DIELECTRIC STRAPS 
SURROUNDING A PERFECTLY CONDUCTING CONVEX 

SURFACE 

2.1 Introduction 

In this chapter, a method is presented to find the electromagnetic fields scat- 
tered by a thin dielectric strap surrounding a perfectly conducting convex surface. 
First, a two-dimensional method is presented to determine the current excited 
across the width of the strap by an incident plane wave. This current is then 
extended to the three-dimensional case using high frequency techniques. Next, 
expressions for the fields radiated by a current element on a perfectly conduct- 
ing convex surface are given using the Uniform Geometrical Theory of Diffraction 
(UTD). Finally, the total scattered field is written in the form of an integral of this 
current over the strap. The integration is then simplified by using a stationary 
phase technique. 

Some of the assumptions which are used throughout this report are the follow- 
ing. The radius of curvature of the surface, S, is large compared to the wavelength 
of the incident wavefront, and the thickness, T, of the strap is much smaller than 
the wavelength. The dielectric material of the strap has a relative permeability, 
Hr — 1, and an ezp(ju>f) time dependence is assumed and suppressed. Whenever 
the words “large” or “small” are used, they mean large or small relative to the 


7 



wavelength. 

2.2 The Current Induced in the Strap 

One way to find the electromagnetic field scattered by a strap is to find the 
current induced within the strap by the incident electromagnetic wave. In this 
section, a method to find that current is described. In this study, a perfectly 
conducting surface, 5, is surrounded by a strap of width, T, and thickness, T, as 
shown in Figure 4(a). 

The first step in evaluating the current in the strap is to divide the strap into 
several smaller sections as shown in Figure 4(b). Then, the current inside each 
section is determined by assuming that the section under consideration is part of 
an infinitely long strap which exists over a ground plane tangent to the surface S 
at the center of the section (see Figure 4(c)). One can see that three assumptions 
are used here. The first one is that each section can be treated alone, which implies 
that there are no interactions between the different sections. The second one is 
that each section exists over an equivalent ground plane. The last one is that the 
current inside a section can be found by considering that section to be part of an 
infinitely long strap. These three assumptions can be justified as follows. 

The first assumption is based on the fact that, diffraction is a local phe- 
nomenon, which means that the scattering from any particular section is simply a 
function of the characteristics of that section only. This implies that the current in 
each section is a function of the characteristics of that section only. Hence, one can 
treat each section alone to find its associated current. The second assumption is 
made because the radius of curvature is large and the sections are small, such that 
each section appears to be mounted over a planar surface instead of a curved one. 
The third assumption is used because each section is part of a continuous strap. 
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DIELECTRIC STRAP PERFECTLY ELECTRIC 



THE DIELECTRIC STRAP DIVIDED 
INTO n SECTIONS 


W 


(b) 



(c) 


Figure 4: (a) The surface, 5, surrounded by the strap, (b) The strap divided into 
sections, (c) Infinitely long strap existing over a ground plane tangent to the 
surface, 5, at the section under consideration. 
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So, the two terminations of each section do not contribute to the current. There- 
fore, the factors which affect the current in the section under consideration are the 
same as the factors which affect the current in an infinitely long strap. Based on 
the Geometrical Theory of Diffraction (GTD), these factors are the incident field, 
the diffraction from the two original edges forming the width of the strap and 
the interactions between these two edges. Therefore, the current in any particular 
section can be found by solving for the current in an infinitely long dielectric strap 
mounted over a ground plane which is tangent to the surface at the center of that 
section. 

Now, let us find the current in the infinitely long strap by using the moment 
method technique. Consider an infinitely long strap mounted on a ground plane, 
see Figure 5(a), and excited by a plane wave which is given by 

j^j _ e fi{X cos 0 sin 6 + Fsin<£sin0 4- ZcosO) ^ 

where X, Y, Z, 0 and <j> are shown in Figure 5, and k is the wavenumber of the 
incident field. The goal is to find the equivalent current inside the strap due to 
the incident plane wave. The first step is to get rid of the ground plane by adding 
the image of the structure. That is, double the thickness of the strap and add the 
image of the incident plane wave. By doing so, one obtains the problem shown in 
Figure 5(b). This problem can be solved by using the method of moments [4,5]. 
In the method of moments, the cross section of the dielectric strap is divided into 
cells which are small enough so that the electric field intensity is nearly uniform in 
each cell. The electric field intensities within each cell are initially considered to be 
an unknown quantities and a system of linear equations is generated by enforcing 
at the center of each cell that the total field must equal the sum of the incident 
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GROUND PLANE 



y 



(c) 


Figure 5: (a) An infinitely long strap over a ground plane, (b) The result of (a) 
after removing the ground plane and adding the image of the incident wave and 
the strap, (c) The strap width divided into cells. 
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and scattered fields. This system of equations is then solved to obtain the electric 
field intensity in each cell. The equivalent electric current in each cell is then given 

by 

^cell n ~ ~ e o) E ce ll n (2*1) 

The above solution requires, as an input, the incident field at the center of 
each cell. For our particular problem the incident field at the strap mounted over 
the surface, 5, shown in Figure 4 is different from the field incident at the strap 
mounted over an infinite ground plane as shown in Figure 5. This difference is 
related to the surface curvature at the location of the strap on the structure. This 
field is found by using high frequency techniques. First, the far field radiation by 
a source located on the surface, S, is found. Then, reciprocity is used to find the 
field on the surface, 5, due to a source located in the far field. This field is used 
to find the current in the short dielectric strap section as described earlier. 

2.3 Scattering from the Strap 

At this point of our development the equivalent electric current in the strap 
is known. The radiated (scattered) electric field can be written using 

E(P)= IJI dE{P/Q') dv' (2.2) 

where dE(P/Q') is the electric field at the observation point, P, due to an in- 
finitesimal electric current moment at the source point, Q 1 . By using UTD [6] one 
can find dE(P/Q'). This is the topic of the next section. 

2.3.1 Radiation From Sources Over Convex Surfaces 

Consider an infinitesimal electric current moment, dp(Q'), located on a per- 
fectly electric conducting convex surface as shown in Figure 6. According to geo- 
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(a) shadow region 


(b) lit region 



metrical optics, the space surrounding the source is divided into an illuminated(lit) 
and shadow regions by a plane tangent to the surface at Q' . This plane is referred 
to as the shadow boundary for that source location. The present formulation of the 
uniform GTD solution leads to separate representations for the elemental radiated 
field dE(P/Q') in the shadow and lit regions, respectively. However, these different 
representations have been shown to match exactly in polarization, amplitude and 
phase at the shadow boundary [7]. 

The unit vectors used in this treatment are shown in Figure 6. The unit 
vector, s, defines the radiation direction. At the source point ( Q '), n 1 is the unit 
surface normal vector, t 1 is the unit tangent vector of the geodesic path, and the 
binormal unit vector is given by V = t' x h'. Similar unit vectors are introduced 
along the geodesic path. Thus at the diffraction point (Q), n is the unit surface 
normal vector, t is the unit tangent vector of the geodesic path and the binormal 
unit vector is given by b = t x n. Note that hi is the unit vector normal to s 
and lies in the plane containing h! and i. Also, i' L is a unit tangent vector to the 
surface, 5, at Q' , and it lies in the plane containing h 1 and i. The binormal unit 
vector is defined by b'i which is given by = t'i x h' or b' L = s x n^. 

In the lit region, one may express the elemental electric field produced by 
dp(Q' ) in a ray optical format as follows [6]: 

dE(P L )~C dp(Q') T L (2.3) 

s 



where — represents the spherical wave associated with the point source at Q 1 . 
The quantity, C, is given by 


C = 


-jk 

47T 


(2.4) 


14 


Likewise, one may also express the field in the shadow region as follows [6]: 

dE(P s )~C dp(Q') • f 


I jgj f EiiSl 1 1 ^ 6 g -jki 

d V (Q)\pg(Q')l 


M A P A 2 A e ~ jksd {2.5) 




The explicit forms of the radiation dyadic transfer functions, T and T, are given 

by [6] 

f L = Z 0 sm6 i {n'h L (H L + r o 2 Fcos0‘) + n'bT 0 F} 


( 2 . 6 ) 


and 


T = Z 0 {n'hH + n'bToS} 


i A 

L r 


(2.7) 


where t is the geodesic arc length from Q' to Q. The distance between the diffrac- 
tion point, Q , and the observation point, P, is given by s d and Z 0 is the free space 
intrinsic impedance. The quantities (H L , H , S L and S ) are referred as the hard 
( H H) and soft (S^; S) Fock functions which describe the radiation patterns of 
a source on a convex surface. These soft and hard Fock functions are plotted in 
Figure 8 and defined by [6]: 

,-ji L r 

5 


r t t°° e~ J<i ' 

H (i >= ^ Le-**I3 dT 


( 2 . 8 ) 


1 too e~ti T 

= ~JS dT ^(7) 


S L (( L ) 




-i( L /3 


ji r dr , and 

It J ooe~3 2 */ 3 Wo(t) 


m 


( Q ') V* 


Sit) — l . 1 f°° 

miQ 1 ) \fit Jooe~i 


dr 


w 2 {r) 


o~j(r 


( Q ') y/ir Jooe-i 2 */ 3 w 2 ( r ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 
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The Fock type Airy function is given by 


W2(T) = 7 ? £-» dT e(T, ~ im 

(2.12) 

with w'^t) being the derivative of u>2( T ) with respect to r. 
£, for the shadow region is given by 

The Fock parameter, 

[ Q jj M (0 

JQ' P g(t') 

(2.13) 

with 


m( f '\ f^ff(^)l 1/3 

m(t ) — g 

(2.14) 


Here Pg(t') is the surface radius of curvature along the ray path at t 1 . The Fock 
parameter, £ L , for the lit region is given by 



£ L — —m L (Q')cosO t 

(2.15) 

with 


m L (Q') = m(Q') [l + T 2 cos 2 0‘] 1/3 

(2.16) 

and cos 0 l = h' 

• i. The torsion factor, T 0 , is given by 



To = T(Q')p„[Q') 

(2.17) 

where T(Q') is 

the surface ray torsion at Q 1 . One may show via 

differential geom- 

etry that [7] 

t(q') — s * n ^ a ( * ) 

2 [ R 2 (Q') Ri(Q') 

(2.18) 

and 

1 cos" a sin* 5 a 

p,m ~ Rim + mQ r ) 

(2.19) 


where R\{Q') and R2{Q') denote the principal surface radii of curvature at Q ' , 
and a' is the angle between t' and f j as shown in Figure 9. Here fj and are the 
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Figure 7: Perspective view of a surface diffracted ray tube (enlarged view), 
principal surface directions at Q' associated with Ri(Q') and f?2 (<?*)> respectively. 


The function, F, is given by 


S L -H l cos0* 

1 + T 2 cos 2 0* ' 


( 2 . 20 ) 


The angles, dip 0 and di/> = are shown in Figure 7. Note that pi j is the wave 

p 2 

front radius of curvature, or the caustic distance, of the surface diffracted ray in 
the b direction at Q and may be calculated via differential geometry. 

Now, let us go back to Equation (2.2) and write it in the following form: 


E(p) = J j ) e l"W ) el''! 1 ,: > kr °d4,'<lt' 


( 2 . 21 ) 


where 


t' is the integration variable along the width of the strap 
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Figure 8: Plots of the radiation Fock functions. 


''i “?i i 

tv t = cos a 



b- 


Figure 9: Principal surface directions at the source location. 
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<f>' is the integration variable along the length of the strap 

e jK$(<l > ) con tains the current and the radiation phase variation along the 
strap 

contains the radiation phase variation along the width of the strap. 

F(<t>' ,t') contains the remaining variables and constants. It does not contain 
any phase variation related to <f>' , and 

r 0 is the distance from the center of the structure to the observation point. 

The triple integral in Equation (2.2) is converted into a double integral in 
Equation (2.21) because the strap is very thin compared to the wavelength such 
that the current is uniform across its thickness. The latter integral can then be 
evaluated numerically; however, the use of the stationary phase method makes it 
much easier to evaluate as described in the next section. 

2.3.2 Stationary phase integration method 

In this section, the stationary phase integration method [8] is summarized. 
Let us consider the following integral: 

/(«) = [ b F(<t>)ei K *Wd<f> ( 2 . 22 ) 

Ja 

where 


<f> is a real variable 
k is real, positive and large 

$(<f>) is a real, continuous function with continuous derivatives for a < <f> < b, 
and 
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F(4>) may be a complex function, but it is slowly varying and well behaved 
in the range of integration. 

Let = 0 for a < <f> a < b, and <f> 9 is not close to a and b. Then, (f> $ is called 

the stationary phase point, and the integral, I, can be approximated by 

/(«) ~ 

+1 _ 1 F ( a ) j[^( a )-f] / 2 23 ) 

«*'(&) «*'(a) { 

where sgn (*) = — 1 if x < 0, and sgn (®) = 1 if x > 0. The above results were 
found by dividing the integral in Equation (2.22) into three integrals and then 
by using an asymptotic approximation to evaluate these integrals. The first term 
in Equation (2.23) corresponds to the contribution from the stationary point, <f> s . 
The other two terms correspond to the contributions from the end points of the 
integration a and b. If the function, $(<£), has more than one stationary point in 
the interval (a, b), one needs to sum the contributions from each of the stationary 
points. If the integration in Equation (2.22) is a closed loop integration then, there 
are no end point contributions. 

In using the stationary phase integration method to simplify the integral in 
Equation (2.21) one needs to split the integral into two integrals. One of these 
two integrals covers the part of the strap which is excited directly by the incident 
wave. The other one covers the part of the strap which is excited by a creeping 
wave in that the strap is shadowed by the structure. Apparently, each integral 
contains one stationary point without an end point. 

In the next chapter, a method to find the scattered electromagnetic field from 
a strap surrounding a perfectly conducting convex structure is developed based on 
the current found in this chapter. 
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CHAPTER III 


BACKSCATTERED FIELDS FROM A DIELECTRIC STRAP 
SURROUNDING A CONVEX STRUCTURE 

3.1 Introduction 

In this chapter, the method described in Chapter II is used to evaluate the 
backscattered fields from a thin dielectric strap wrapped around a structure such as 
shown in Figure 10. Recall that the relative permeability of the strap is assumed 
to be unity; while, the dielectric constant of the strap is e r . For simplicity the 
pointing vector of the incident wavefront lies in a principal plane of the structure; 
i.e., the XZ plane in Figure 10. The backscattered fields are computed in the E- 
plane (XZ plane in Figure 10). The H-plane backscattered fields from a thin strap 
wrapped around a structure are extremely low and one does not need to worry 
about their effect on the a measurement in that the worst case E-plane scattered 
fields set the measurement error level. This claim is justified at the end of this 
chapter. 

The techniques developed here can be used to treat an arbitrary convex struc- 
ture; however, the remaining discussion will be devoted to the ogive. The ogive 
is interesting in that it represents a low level scattering target; thus, one can ob- 
serve the effect of the strap scattering for a practical case. This will become more 
apparent in the next chapter, when the results are examined. 

In the case of a strap wrapped around an ogive there are two different mecha- 
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nisms leading to the currents excited on the strap and then to the scattered fields 
from the current elements. The first one is the direct excitation which occurs in the 
lit part of the strap. The second one is the creeping wave excitation which occurs 
in the shadowed part of the strap. These two mechanisms are treated separately. 

3.2 Description of the Ogive 


The ogive is used in this study to represent a general convex structure. It is 
shown in Figure 10 and can be described by knowing the two dimensions of the 
ogive, a and 6, as shown in Figure 11. The vector R((j),Z) which connects the 
origin and a point on the surface of the ogive is then given by 


Z) = A(Z) cos (<f>)z + A{Z) sin {<f,)y + Zz (3.1) 


where 

A(Z) = -c+sjr 2 - Z 2 

b 2 + a 2 

r = — , and 

2a 

c — r — a . 

Another, convenient, way to describe the ogive is the following: 

R{<t > , v) = (r cos(i^) — c) cos(<^)x + (r cos(t') — c) sin(<^>)$ + r sin(^)z 


(3.2) 

(3.3) 

(3.4) 


(3.5) 


where u is as shown in Figure 11, and from which one can see that 


sin(t') = — . 


(3.6) 


The advantage of the second definition is that the derivative of R with respect 
to v is tangent to the surface of the ogive and perpendicular to <j>. Note that <j> 
is also tangent to the surface of the ogive at any point. So, the two tangents 
perpendicular to each other at any point on the surface of the ogive can be easily 
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found. These two tangents are enough to define the tangent to the surface at any 
direction. They are as follows: 



3.3 Backscattering From the Strap 

Following the method described in Chapter II let us set up the integral for 
the radiation by the current induced in the strap as follows: 

E(p) = [ 2 * [ L/2 F(<f>\t')ej K *W'')ej t 'Pe-j kr °d<f>'dt’ . (3.8) 

J 0 J — L f 2 

The different constants and variables which appear in the above equation were 
defined earlier in Chapter II. The stationary phase integration method is used to 
evaluate the above integral with respect to <f> f . The stationary points are located 
at the top and the bottom of the ogive; i.e. the points closest to and farthest 
from the transmitter/receiver antenna. The one at the bottom is shadowed by the 
ogive; while, the other one is not. This makes the integral formulation at each 
point different from the other one. Therefore, one must deal with each stationary 
point using a different method. 
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3.3.1 The Lit Part of the Strap 


In order to find the backscattered field from the strap, one needs to find the 
current within the strap excited by the incident plane wave. In order to find that 
current let us use the moment method approach described in Section 2.2. That 
method requires, as an input, the value of the incident field at the location of each 
cell in the strap. This can be done as follows. Consider a current moment located 
at the top of the ogive and radiating in free space as shown in Figure 12. The 
radiated field can be found using the formulas given in Section 2.3.1. So for an 
observation point located in the X Z plane, the radiated field is given by 

_ - r p-jks 

dE ~ -r~dp{Q) • f (3.9) 

47T s 

where 

f L = Z o sin0 i H L {-7n(Q')cos0 i )h'h L , and (3.10) 

MQ') = dp(Q')n' . (3.11) 

Then, one finds that 

dE~ ^Z o dp(Q , )sin0 i H L (-m(Q')cos0 i )^-^-h L . (3.12) 

47 r s 

The different variables in the above equation were defined in Chapter II. Now, 
by using reciprocity, one can find the field on the surface due to a plane wave 
incidence. By doing so, the electric field at the surface, 5, is given by 

H L (-m(Q') cos0 i )sm0 i e^ kasin9cos< f , 'e^ Zcos0 ) . (3.13) 

The above equation is used in the moment method approach to find the current 
in the strap. That current is then used to find the backscattered field from the lit 
part of the strap. 
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The radiated electric field by the current in the lit part of the strap is as 
follows: 

- ff J-Z I-l • r L€ ^° d +' di < 314 > 

where 

Z Q is the intrinsic impedance of free space 
T is the thickness of the strap 
k is the wave number of free space 

R is the distance from the source point to the observation point 

<f>b is the angle which defines the separation boundary between the shadowed 
and the lit part of the strap, and 

J(<f>',t') is the current in the strap. 

Equation (3.14) can be written as 

-ik e~i kro fL/2 r<t> h T 

E U (P) = — Z 0 T / [ sin O t H L (—m(Q) cos 0 l ) 

4ir r 0 J-L/2J-<f>i, 

J(<f>' ,t' )h L e^ a sin ^ cos ^ ) e ^ kt ' cos ® ad^'dt 1 (3.15) 

where 

6 is the observation angle measured from the Z— axis, and 
r 0 is the distance from the ogive to the observation point. 

Now, the current can be described by 

= J 1 (<f> , ,t l )e? ka,inecoB *' (3.16) 
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Then, Equation (3.15) can be written in the following form: 


/•i/2 

E u (p) = / 

J-L/2 4 


(3.17) 

F(<f> , ,t')=^-Z 0 sin 
4nr 0 

ie i H L (-m(Q') cos O^TaJ^’ ,i')h L 

(3.18a) 


k = 2fcasin0 

(3.186) 


$(^) = cos <j> , and 

(3.18c) 


f3 = k cos 6 . 

(3.18d) 


Now, the above integral can be evaluated using the stationary phase integra- 
tion method. First, the stationary phase point is found from 


4>{<i>) = cos 4> 


(3.19a) 


$ , (^ ; ) = — sin^ , and 


(3.196) 


$"(<£') = -cos<j) (3.19c) 

where the stationary point is given by 

*W) = 0 • (3.19d) 
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In this case, the stationary points are given by <j> a t = 0,7r. Since the interest here 
is in the lit part of the strap, <f> a i = 0 is the only stationary point to be evaluated 
here. Using the results of section 2.3.2, one finds that 


E 


V) = f_ 


L/2 _ 


2ir 


I 


e j(^(4>3tHisgn(^"(4> st )) e jt'(3 e -jk r o dt i (3 20 ) 


which can be written as 

E u (p) = e~ ikr ° f L/2 F(0,t')e ji,(3 dt' . (3.21) 

V ka sin 0 J-L/2 v ’ 

The above integral can then be easily evaluated numerically across the width of 

the strap. 

3.3.2 The Shadowed Part of the Strap 


According to Equation (2.5) the radiation by a current element in the shad- 
owed part of the strap can be written as 


E S {P) = 


- jkT 


rh/i ti 

4tt J-L/2 J<t> b 


dip 0 

PgiQ) 

1/6 •, 

e~ jkt 

P 

.Pg(Q'). 

N 


P 2 


N Api + > d ) 


; * ka ad<t> dt 1 


(3.22) 


where 


J (</>', t') is the current in the strap. This current is found in a similar way 
to that used for the current in the lit part of the strap. The details are 
given at the end of this section. The other factors in Equation (3.22) 
are defined as follows: 


T = Z Q h'hH(£) + n , 6T 0 5(^)j , but because of the symmetry the b compo- 


nent cancels out; therefore, one needs to consider 


f = Z o h'nH{0 


(3.23a) 
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and p 2 are the ray divergence factor and the caustic distance, respec- 
tively. They are found using the cylinder perturbation method [9]. See 
Appendix A for details. 

Pg(Q) and Pg{Q') are the surface radii of curvature in the t and t' direction, 
respectively. They are given by Equation (2.17). 

t is the geodesic path length, which is given by 

t = — rO cos <f> — t 1 (3.23 b) 

t 1 and 4> define the position of the source, which are shown in Figure 13. 

s d is the distance between the diffraction and the observation point, which 
is given by 

s d = r a + (r — a) sin 9 cos <f> (3.23c) 

r 0 is the distance from the center of the ogive to the observation point, and 
finally 

£ is given by 


t = f Q 

JQ’ Pg(t') 


, and 


(3.23d) 


m ( Q . )= p £^] I/3 


(3.23e) 


Using this approach, the shadow region field can be written as 

E\P)=f L/ 2 f 2 ^ F((A^/')e• ,K# ^' ) e^c-• ?7^ro d<A , df , (3.24) 

J—L/2 J<bu 


where 




\ P2 d i' \ sd \p2 + sd ) 


(3.25a) 
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wave radiation from a source mounted on an ogive. 
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k = 2 (r0k — k(r — a) sin 0) 


(3.25 6) 


$(< f /) = cos <f> (3.25c) 

f3 = k , and (3.25 d) 

= /(</,', <')e-J (r ^- fc ( r ^ a)sin0) . (3.25e) 

The function H(£) has a phase dependence on <f>'] however, it is small, and it will 

not affect the stationary phase integration result. 

Now, let us use the stationary phase integration method to evaluate the in- 
tegral in Equation (3.24) with respect to <f >' . First, one finds the stationary phase 
point using the following expressions: 

$(<^/) = cos <f> (3.26a) 


$\<}>) = — sin<^ , and (3.266) 

and 

$"(<£') = -cos <£' . (3.26c) 


The stationary point is at ^{4> s t ) = 0 which means that <f> 3 i = 0,ir. But since the 
shadowed part of the strap is being studied here, <f> 3 t = ?r. Then, oiie finds that. 


rL/2 _ 

' iW) 1 


k|$ m (7t) 


J[«M "•)+ 1 sgn($"(jr))] e jt'P e -jkr odi i (3 27) 


33 


which can be written as 


E a (p) = y^ e J (- K+ ?) e --'* r ‘> J L/ L 2 ^F(ir,t')e jt 'Pdt' . (3.28) 

The above integral is then evaluated numerically. 

The current is found using the moment method approach described 

in Section 2.2. In using that method, one needs to find the excitation for each 
cell in the strap; this can be done as follows. Consider an electric current moment 
located at the bottom of the ogive radiating in free space as shown in Figure 12. 
According to Equations (2.5), (2.19) and (3.23a), the electric field due to a current 
moment at the stationary point and an observation point located in the X Z plane 
is given by 

(3.29) 

(3.30) 
(3.316) 

s* = r 0 — (r — a) sin 6 . (3.31c) 


dE(p a ) ~ C dp(Q')Z o H(0^ ^ d~i ki t 


, p 7 e~^ d n 

AP2 + 3 ) 


In the far field, the above equation can be written as 


dE(p a ) ~ C dp(Q )Z 0 H(£)J -~e 


tyo-jkt 1 r -7 ks d ~ 


dr!) $ C 


t = rd — t f , and 


One can use reciprocity to find the electric field at the location of the electric 
current moment due to a plane wave incidence. It is found to be given by 



— jkt 1 gj{ — rQk+k(r—a) sin 6) 


(3.32) 
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The total backscattered electric field is then given by 

E(P) = E U (P) + E S {P) . (3.33) 

In Equations (3.21) and (3.28), one can see that, when 0 — 0°, the field 
becomes singular which means that the equation is not valid at or near that caustic 
direction . This is true because at 9 — 0°, all the current elements around the 
strap contribute to the backscattered field. In order to find the field at or near this 
caustic direction, one needs to integrate the current around the whole strap. This 
can be done as follows. When 0 = 0° the electric field radiated by the current in 
the strap can be written as 

-ik t r L / 2 , , e~i kR , . 

B(r) = -~ZoH L (0)T j_ J(4/,t')h L —^-ai4Ut! (3.34) 

Let R = r 0 — i’\ then, one finds that 

E(p) = ^-Z 0 TaH L { 0)e~i kr ° J(4>', t')n L e jkt ' d<f>'dt' . (3.35) 

From the symmetry of the structure one can see that the backscattered electric 
field in * direction can be written as 

E(p) = ^-Z o TaH L {0)e- jkr ° \ L> 2 [** J(<t>\t')cos(<f>')e jki ' d^di'x (3.36) 

47 tVq J — Jj f 2 JO 

Also, one can easily see that J (<(>', t 1 ) can be written as 

J(<f>',t') = J(0,t')cos(<f>') (3.37) 

Thus, one obtains that 

E{p) = ^-Z 0 TaH L (Q)e~ jkr ° [** J{0,t') cos 2 (<f>' )e^ t ' d^'dt'x (3.38) 

47rf(j J — i/2 Jo 

and 

E(p) = ^Z o TaH L (0)e-J kr ° j L ^J(0j')e jkt ' dt'x . (3.39) 
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The above integral can again be evaluated numerically across the width of the 
strap. 

Let us now justify the claim stated at the beginning, the backscattered H- 
plane pattern is very low and does not need to be considered. The backscattered 
field from the strap described in equations (3.21) and (3.28) depends on the elec- 
tric current at the stationary points. That current is a function of the incident 
electric field at the stationary points. In the H-plane backscattered field case, the 
stationary points are at <j>' = tt/ 2 and <f>' — —ir/2. The electric field is very weak at 
these stationary points because the incident electric field is tangent to the surface 
of the, and the strap is very thin. This means that the electric current at <f> — tt/2 
and <f>' = —tt/2 is almost zero, hence the backscattered field in the H-plane is very 
low. 

In this chapter, the method presented in Chapter II was used to solve for 
the backscattered fields from a strap surrounding a structure. This technique was 
applied to an ogive simply to provide results for a specific case. The method split 
the solution into two parts. One part corresponds to the backscattered field from 
the lit part of the strap, and the other one corresponds to backscattered field 
from the shadowed part of the strap. The above method does not work at or near 
6 = 0°, because it is a caustic for the backscattered field such that stationary phase 
dose not apply. Therefore, at that point the backscattered field was evaluated by 
integrating the current over the whole strap. In the next chapter some computed 
and measured results are presented. 
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CHAPTER IV 


RESULTS AND EMPIRICAL DESIGN FORMULAS 

4.1 Measured and computed results 

In this section, the results of chapter III are used to compute the E-plane RCS 
of a strap wrapped around an ogive as shown in Figure 14. First, some computed 
and measured results for a particular strap are compared in order to demonstrate 
the validity of this theoretical approach. Then, some individual computed terms 
are presented to demonstrate the way in which the solution works. 

Figure 15 shows the RCS of a dielectric strap wrapped around an ogive versus 
frequency for nose-on incidence (0 — 0°). The strap has L — 2 ", T = (hOtM” and 
e r = 2.3. The dimensions of the ogive are a = 2.37 M and b = 18 ,f . Note that 
at 4GHz the RCS of the strap is about -65 dBsm. The RCS increases with the 
increase in frequency; however, it is still low at high frequencies. In Figure 15, the 
measured RCS is also shown which show reasonable agreement with the calculated 
ones. To measure the RCS of the strap, it was wrapped around an ogive and 
mounted in the OSU compact range facility as shown in Figure 16. First, the 
RCS of the whole structure was measured from 2 to 18GHz every 10MHz. Next 
the phase center was taken to be at the middle of the strap; then smoothing was 
used to isolate the RCS of the strap. The measured result which appears in Figure 
15 was obtained using a 201 point smoothing. A 101 points smoothing gives the 
result shown in figure 17. One can see that the 201 point smoothing gives the most 
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Figure 15: Nose-on RCS of a strap versus frequency with 201 point smoothing. 
a=2.37", b=18.0", L=2", T=0.004" and 6 r =2.3. 
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Figure 17: Nose-on RCS of a strap versus frequency with 101 point smoothing. 
a=2.37", b=18.0", L=2", T=0.004" and c r =2.3. 
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reasonable results. Therefore, the rest of the measured results are presented using 
a 201 point smoothing. The time domain results obtained from the measured data 
is shown in Figure 15(a). It indicates that the strap scattering is dominated by 
the front and back edges. Note that the two edges of the strap scatter roughly the 
same level, which causes the lobe structure observed in the frequency domain. 

Figure 18 shows the measured as well as the computed RCS of the tip of the 
ogive. The computed RCS of the tip of the ogive was found by using the following 
physical optics formula: 

-(°)=^2 tan V) (4- 1 ) 

where a is the half cone angle of the ogive. One can note that the measured 
RCS does not agree with the computed one near the upper and lower ends of the 
frequency spectrum because the measured result was obtained using the smoothing 
approach mentioned earlier. Note that smoothing does not work well at the ends 
of the spectrum because one can’t center the smooth data around the ends. In 
this case, the end values are found using 101 points to the right or left of the end 
point. 

Figures 19 through 26 show the computed as well as measured RCS of the strap 
versus frequency for observation angles from 5 to 40 degrees with an increment of 
5 degrees. The quality of the measured results decreases with the increase in the 
observation angle, 9 , in that as 9 increases the return from the strap gets smaller 
and smaller. This make it hard to separate the RCS of the strap from the noise. For 
Figures 19 through 26, the RCS was computed without including the contribution 
from the shadowed part of the strap in that the cylinder perturbation method does 
not work well for thin ogives, such as the one used for the measurement. However, 
this did not have any considerable effect on the results. This point well be discussed 
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Figure 18: Nose-on RCS of the tip of the ogive versus frequency with a=2.37" 

and b=18.0". 

later in this section. Although the computed results in Figures 19 through 26 were 
found using this simple approach, reasonable agreement was obtained in each case. 

In order to demonstrate how the solution works, some computed results for 
different straps are presented. Figure 28 shows the backscattered field from a 
dielectric strap over a ground plane, with the geometry shown in Figure 27. The 
backscattered field is computed at 6GHz, and the dimensions of the strap are 
L = 2.5", T = 0.004" and = 4.0. The null which appears at 9 = 90° is expected 
for a thin strap because the incident electric field and its image field cancel at the 
ground plane, which means that the strap is not significantly excited and hence 
does not backscatter significantly. A similar null is expected in the RCS pattern 
of the strap wrapped around the ogive. 

Figures 29 through 31 are for an ogive with dimensions, a = 14.76" and 
b — 29.53". Figure 29 shows the RCS of a strap with L = 3.69", T = 0.008" and 
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Figure 19: RCS of a strap versus frequency with a=2.37”, b=18.0", L=2”, 

T=0.004" and e r =2.3. 



FREQUENCY IN GHZ 


Figure 20: RCS of a strap versus frequency with a=2.37 ,/ , b=18.0 ,/ , L=2 M , 

T=0.004" and e r =2.3. 
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Figure 21: RCS of a strap versus frequency with a=2.37 /; , b=18.0 ,, ) L=2 

T=0.004" and e r =2.3. 


I 

CD | 
Q ' 
o 

LJ o 

si 


8.0 10.0 12.0 14.0 16.Cfc> 


RCS ENVELOPE 


# = 20 ° 



14.0 


:ii > -g 

MEASURED \j T 

CALCULATED U p 

6.0 8.0 10.0 12.0 14.0 16.01 

FREQUENCY IN GHZ 


Figure 22: RCS of a strap versus frequency with a=2.37 ,, > b=18.0 ,, ) L=2 

T=0.004" and e' r =2.3. 
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Figure 23: RCS of a strap versus frequency with a=2.37", b=18.0", L-2", 

T=0.004" and e r =2.3. 
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Figure 24: RCS of a strap versus frequency with a=2.37", b=18.0", L-2" 

T=0.004" and e r =2.3. 
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Figure 25: RCS of a strap versus frequency with a=2.37 ,/ , b=18.0 M , L=2 ,/ , 

T=0.004" and e r =2.3. 
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Figure 26: RCS of a strap versus frequency with a=2.37 ,/ , b^lS-O^, L=2 ,, ) 

T=0.004" and e r =2.3. 
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0.004” and e r =4.0 and 2.0 . Figure 29(a) shows the RCS of the lit part of the 
strap; whereas, Figure 29(b) shows the RCS of the shadowed part of the strap. 
The total RCS of the strap is shown in Figure 29(c). Figure 30 shows the RCS of 
a strap with L = 2.36 ,, J T = 0.004” and 0.008”, and e r = 4.0. As mentioned in 
chapter III, one can clearly see that the RCS of the lit and shadowed parts of the 
strap are not valid at or near 9 = 0°. Therefore, at 6 — 0° the RCS of the strap is 
evaluated by integrating the current all around the strap. In the total RCS of the 
strap the, results are only shown outside the caustic region. 

Figures 29 and 30 show that the RCS of the strap can be reduced by reducing 
the thickness and/or the dielectric constant of the strap. This means that a tradeoff 
can be made between the RCS of the strap and its physical properties, such as 
the maximum load that it can handle which is influenced by the thickness of the 
strap and its dielectric constant. An empirical formula which can help in making 
such a tradeoff is given in [2]. One can easily see that the nulls which appear in 
Figures 29 and 30 (except for the one at 0 = 90°) are due to the contributions 
from the front and back edges of the strap coming out of phase. Figures 29 and 30 
also show that changing the width of the strap does not affect the RCS level but 
it affects the number of nulls appearing in the RCS pattern. 

Figure 31 shows the nose-on RCS of the strap versus frequency. The RCS 
level of this strap is higher than the one in Figure 15, the main reason being that 
the dimension ‘a’ in the ogive of Figure 31 is equal to 14.76”; while, ‘a’ in Figure 
15 is equal to 2.37”. This means that the strap of Figure 31 is much longer than 
the one of Figure 15, which explains the difference in the RCS level. 

Figure 32 shows the RCS of the same strap used in the measurement (i.e. 
L = 2”, T = 0.004” e r = 2.3, a = 2.37” and b = 18”) at a frequency of 4GHz. 
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(c) Total RCS of the strap. 



Figure 29: Continued. 
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(c) Total RCS of the strap. 


Figure 30: Continued. 
















The creeping wave contribution to the RCS of the strap was not included because 
the cylinder perturbation method gives inaccurate results for such a thin ogive in 
that for a thin ogive the increase in the length of the creeping wave path with the 
increase in the observation angle, 9 is larger than what the cylinder perturbation 
method can handle. Therefore for thin ogives, the cylinder perturbation method 
gives good results only for small 9. But for small 9, the solution for the shadowed 
part of the strap is not valid because of the caustic mentioned earlier. Another 
perturbation method such as the cone perturbation method [9] could be used. 
However, one will not gain too much from that, because the creeping wave term 
will stop at 9 = 15° when it hits the tip of the ogive. Therefore, the RCS of the 
strap is computed at 9 = 0° then the RCS of the lit part of the strap is computed 
alone beginning at 9 = 2 or 3 degrees (i.e. outside the caustic region). One should 
not pay much attention to the discontinuity appearing in the RCS pattern, because 
one is interested in the RCS envelope rather than the exact RCS pattern. Figure 
33 shows the RCS of the same strap at 10GHz. The envelope of the RCS of the 
strap increased by about 8 dB with the increase of the frequency from 4 to 10 
GHz. This can also be seen in the frequency domain plots as shown in Figure 15. 

4.2 Empirical Design Formulas 

In this section, the computed results given in the previous section and the 
analytical expressions given in chapter III are used to establish some formulas 
which relates the RCS of the strap shown in Figure 34 to the different parameters 
of the problem. Such formulas are useful in designing and using the straps as a 
target support structure. 

Lai and Burnside [2] provided a formula which relates the RCS of a strap in 
free space to its thickness, dielectric constant and the operating frequency. That 
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formula stated that the RCS of the strap has (T) 1 " 7 x (e r ) 103 dependence. It can 
be seen that the results given in the previous section has the same dependence. 

For 0 S — 0, Equation (3.39) along with the computed results of Section 4.1 
are used to obtain the following equation: 


<r 5 (0) = 0.833 /- 53.084 + 10 log(e r ) 103 + 10 log(T) 1 - 7 

+ 20 log(L,) (4.2) 


where 


<r a (0) is the RCS of the strap in dBsm at 6 S = 0, and 0 S is the incidence 
angle as shown in Figure 34 

/ is the frequency in GHz 

T is the thickness in inches 

L a is the strap length, for example L s = 27ra for the strap in Figure 34, and 
e r is the dielectric constant. 

For 0 $ outside the caustic region, say 6 S > 3°, Equation (3.21) along with the 
results of Section 4.1 are used to obtain that 


<r 9 (d 3 ) = 0.5/- 49 + 10 log(c r ) 103 + 10 log(T) 1 ' 7 

+ 20 log(cos(0 4 )) + 10 log( R ) a ) + 401og(^ I (O) (4.3) 

2s si | 

with 

( = (~ ~^-- )1/3 sin(^) (4.4) 

where 
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Qs is the stationary point in the lit part of the strap, it is taken to be in 
the middle between the two edges of the strap 

R i is the radius of curvature of the surface of the structure at Q a in a 
direction parallel to the edges of the strap 

i?2 is the radius of curvature of the surface at Q s in a direction perpendicular 
to the edges of the strap. 

H^(£) is the Fock radiation function, it is defined in chapter II, and 

9 a is the angle between the observation direction and the principal plane of 
the strap, as shown in Figures 34 through 36. 

Note that the RCS envelope of the strap does not depend on the width of the 
strap; therefore, the width does not appear in Equations (4.2) and (4.3). Also, 
note that Equations (4.2) and (4.3) are good for general convex structures, and 
the ogive is chosen for demonstration purposes only. 

Now, using Equations (4.2) and (4.3) the RCS envelope of the strap versus /, 
9 a , T , e r , JJj, and i?2 ma y t> e found. Equation (4.2) is used to find the RCS at 
9 a =0°, and Equation (4.3) is used to plot the RCS envelope for 9 a outside the 
caustic region, say 9 a > 3°. For 0° < 9 a < 3°, the RCS plot can be simply taken 
to be a constant value. 

Figures 37 through 39 show the RCS envelope of the strap shown in Figure 
34 with a = 2.37" and b = 18 w . Figure 37 shows the nose-on RCS envelope of the 
strap versus frequency with T = 0.002”, 0.004” and 0.006”, and e r = 2.3. Figure 
38 shows the RCS envelope with e r = 2, 3 and 4, and T = 0.004”. Figure 39 shows 
the RCS envelope of the strap versus the observation angle, 9 S with T = 0.002” 
and 0.004”, and e r = 2.0 and 4.0. 
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Figure 35: Ogive hung by straps. 
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Figure 37: Nose-on RCS envelope of a strap wrapped around an ogive versus 
frequency with a = 2.37", b = 18", = 2.3 and T = 0.002", 0.004" and 0.006". 
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Figure 38: Nose-on RCS envelope of a strap wrapped around an ogive versus 
frequency with a = 2.37", b = 18", T = 0.004" and e r = 2, 3 and 4. 
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Figure 39: RCS envelope of a strap wrapped around an ogive with a = 2.37”, 
b = 18", T = 0.002", 0.004" and e r = 2, 4 and f=4GHz. 
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Now , Equation (4.3) is used to estimate the the RCS envelope of the straps 
shown in Figure 35. It is seen that the straps are oriented in a way which minimizes 
its RCS. That orientation makes the backscattering field from the straps appear 
to be coming from two stationary points, these two stationary points are shown in 
Figure 35. One can see that the dominating contribution to the RCS of the straps 
is coming from the leading stationary point because it lies in the lit region; while, 
the other stationary point lies in the shadow region. 

Figure 40 shows the RCS of the leading strap with a = 5 " , b — 20” and 
Z a = 10”, where Z 9 defines the stationary point location. Using the equations 
given in Section 3.2, one finds that R\ = 3.8”, i?2 = 42.5” and 6 S = 13.6°. The 
RCS of the leading strap is plotted for T = 0.002” and 0.004”, and e r = 2.0 
and 4.0. The contribution from the stationary point in the shadow region is not 
expected to have significant effect on the total RCS; therefore, the plot in Figure 
40 is considered to be a reasonable first-order estimate for the RCS envelope of the 
straps shown in Figure 35. 

The basic conclusions of this report are given in the next chapter. 
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Figure 40: RCS envelope of the strap shown in Figure 35 versus frequency with 
a = 5", b = 20", Z, = 10", T = 0.002', 0.004" and e r = 2,4. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 


A method to calculate the electromagnetic scattered fields from a thin dielec- 
tric strap wrapped around a convex, perfectly conducting surface was presented. 
This was accomplished using a moment method technique to find the current ex- 
cited within the strap by an incident plane wave. Then, Uniform Geometrical 
Theory of Diffraction, along with a stationary phase integration method was used 
to compute the fields scattered by the strap. This method was used to evaluate 
the RCS of a strap wrapped around a convex conducting structure. The computed 
RCS values showed reasonable agreement with measured results. It was shown that 
for low frequencies, say below 1GHz, the strap scattering is very low. It increases 
with an increase in frequency, but it is still low at high frequencies. The RCS of 
the strap can be reduced by reducing the thickness or the dielectric constant of the 
strap. So a trade off can be made between the RCS of the strap and its physical 
properties. The H-plane RCS pattern is known to be much lower than the E-plane 
pattern; therefore, only the worst case, the E-plane pattern, was presented in this 
study. 

Some empirical design formulas, which can help in designing the straps, are 
also given. These formulas are good for general convex structures and are expected 
to give resonable estimates for the RCS of the dielectric straps when used as a 
target support structure. 
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APPENDIX A 


RAY DIVERGENCE FACTOR AND CAUSTIC DISTANCE 


The ray divergence factor used in Equation (3.22) can be found as 

follows. First, we define the source location and the radiation direction. Next, we 
find the geodesic path that the ray will travel along it to reach the diffraction point, 
let the path length be L. At the diffraction point the radiation direction i must 
be tangent to the surface and aligned with the geodesic path as shown in Figure 
41. The next step is to draw another geodesic path which has a length equal to L 
and passes through the source location making an angle dj with the first geodesic 
path. The last step is to find the vector i a which is tangent to the surface and 
aligned with the second geodesic path at its end. Then the ray divergence factor 
is given by 


iMQ‘) 


d 7 


(A.l) 


\ diJ>(Q) ^cos -1 (£-< 0 ) 

The above steps can be implemented numerically [10]. However, another 
approximate, but easier, method is followed here. This method is called the per- 
turbation method [9]. The cylinder perturbation method is used here. It is good 
enough when the source exists near the region Z = 0 on the ogive. The method 
uses some of the equations governing the geodesic path on a cylinder assuming 
that they are valid for the ogive. 

The circular cylinder shown in Figure 42(a) is a developed surface; therefore, 
the geodesic path, Q'Q , is a straight line on the unfolded planer surface. As shown 
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GEODESIC PATH 


Figure 41: Creeping wave radiation from source over an ogive 
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(b) Unfolded Planar Surface 


Figure 42: Geodesic path on a developed cylinder. 
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Figure 43: Circular cylinder perturbation, 
in Figure 42(a) the geodesic unit tangent t is given by 

t = cos 7 <f> + sin 7 z (A. 2) 

where 

<}> = ~ sin <f> x + cos <j> y . (A. 3) 

An important thing to put in mind is that, along a given geodesic path the angle 
7, shown in Figure 42(b), is constant. 

Now, as shown in Figure 43, the cylinder is perturbed by bending it around 
the ogive and the geodesics are to be resolved on this perturbed cylinder model. 
Then, using the fact that 7 is constant along a given geodesic path, Q'Q, on the 
perturbed cylinder one obtains the geodesic equation 

tan 7 = — (-4*4) 
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Figure 44: The different parameters for the geodesic path, 
where s e and s r are as shown in Figure (A.4). And 


Sj* — a<j> . 

The point Q is defined by v and <f>. v is equal to the observation angle 0, and <f> is 
taken to be zero. Now, since the cylinder is a developed surface 


A 2 - 4 


(A. 5) 


Then 


cot 7 = 


a<f> 


(A.6) 


^/.L 2 — a 2 <^ 2 

By differentiating both sides of the above equation with respect to <j> and after 
simplification one gets 

(xL ^ 1 

— csc 2 {'f)d') = d<t> . (A. 7) 

s e 
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d<j> can be found from the above expression for any dy. 


Now, one can simply see that 


let 


s e = ru — t' 


jJ- 

r 


then, one gets 

s e = r ( u - u') 

To find du due to a small dy we use Equation (A. 4). 

s e r(v - u') 
tan 7 = — = — — - — - 

S f> Cl(j) 

tan ( 7 + it) = + dt) ] > and 

du = -d<f>ta.n(y + dy) — u + u' . 


(A.8) 

(A.9) 

(A.10) 

(A.ll) 

(A.12) 

(A.13) 


Now, the adjacent diffraction point Q a : (<f> + d<f>, u + du) is known and 


i a = f® cos(7 + dy) 4- f® sin(7 + dy) 


(A.14) 


t — t r cos(7) + t e sin(7) (A. 15) 

where t and t e are functions of the point location, they are given in Section 3.2. 
The ray divergence factor is then given by 


1 


dj>c,(Q') 

dipiQ) ^ 


dy 


cos *(t • ia) 


(.4.16) 


where the angle dy can be conveniently chosen as 1°. 

The caustic distance is easily found. It is the distance between the diffrac- 
tion point ( <f > , u) on the ogive and the intersection point of t and t a when they pass 
through (4>,u) and (<f> + d(f>,u + du ), respectively. 
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